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I Abstract 
-I— > ' 

C ' In July 1925 Heisenberg published a paper [Z. Phys. 33 879-893 

(1925)] which ended the period of 'the Old Quantum Theory' and ush- 
ered in the new era of Quantum Mechanics. This epoch-making paper 
is generally regarded as being difficult to follow, perhaps partly because 
Heisenberg provided few clues as to how he arrived at the results which 
he reported. Here we give details of calculations of the type which, we 
r> | suggest, Heisenberg may have performed. We take as a specific example 

^ . one of the anharmonic oscillator problems considered by Heisenberg, and 

use our reconstruction of his approach to solve it up to second order in 
perturbation theory. We emphasize that the results are precisely those ob- 
tained in standard quantum mechanics, and suggest that some discussion 
of the approach - based on the direct computation of transition ampli- 
tudes - could usefully be included in undergraduate courses on quantum 
mechanics. 



1 Introduction 

Heisenberg's paper of July 1925 1 , entitled 'Quantum-mechanical re- interpretation 
of kinematic and mechanical relations' 2,3 , was the breakthrough which quickly 
led to the first complete formulations of quantum mechanics 4 ' 5,6 . Despite its 
undoubtedly crucial historical role, Heisenberg's paper is not generally referred 
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to in undergraduate courses on quantum mechanics - in contrast, say, to the 
place of Einstein's 1905 paper in the teaching of Relativity. Indeed it is still 
widely regarded as being difficult to understand and - perhaps because of this - 
of only historical interest today. For example, Weinberg 7 has written as follows: 
'If the reader is mystified at what Heisenberg was doing, he or she is not 
alone. I have tried several times to read the paper that Heisenberg wrote on 
returning from Heligoland, and, although I think I understand quantum me- 
chanics, I have never understood Heisenberg's motivations for the mathemati- 
cal steps in his paper. Theoretical physicists in their most successful work tend 
to play one of two roles: they are either sages or magicians.... It is usually not 
difficult to understand the papers of sage-physicists, but the papers of magician- 
physicists are often incomprehensible. In that sense, Heisenberg's 1925 paper 
was pure magic. 

Perhaps we should not look too closely at Heisenberg's first paper ' 

There have, in fact, been many discussions aimed at elucidating the main 
ideas in Heisenberg's paper, of which Refs 3, 8, 9-18 no doubt represent only 
a partial selection; of these the most detailed appear to be Ref. 3 pages 28- 
35, Ref. 8 pages 204-224, Ref. 10 pages 161-188, and Ref. 11 chapter IV. Of 
course, it may well not be possible ever to render completely comprehensible the 
mysterious processes whereby magician-physicists 'jump over all intermediate 
steps to a new insight about nature' 19 . In our opinion, however, one of the 
main barriers to understanding Heisenberg's paper, for most people, is a more 
prosaic one: namely, that he gives remarkably few details of the calculations 
he actually performed, in order to arrive at his results for the one-dimensional 
model systems which he treats (anharmonic oscillators and the rigid rotator). 

One aim of the present paper is therefore to make Heisenberg's paper more 
accessible to scientists, and to historians of science, by briefly reviewing the 
line of reasoning he followed in setting up his new calculational scheme (section 
2), and then by presenting (in section 3) full details of a calculation typical 
of those we conjecture that Heisenberg himself performed. Our 'reconstruc- 
tion' is based on the assumption that, having formulated a scheme which was 
capable in principle of determining uniquely the relevant physical quantities 
(transition frequencies and amplitudes), Heisenberg then applied it rigorously 
to various 'toy' mechanical systems, without any further recourse to the kind 
of 'inspired guesswork' that characterised the Old Quantum Theory. Surpris- 
ingly, this point of view appears to be novel: MacKinnon 10 , and Mehra and 
Rechenberg 11 , for example, consider that Heisenberg arrived at the crucial re- 
cursion relations (equations (|4*U)l - (|4^|) below), in the quantum case, by essen- 
tially guessing the appropriate generalisation of their classical counterparts (see 
section 3). We are unaware of any evidence, now, that can settle the issue. In 
any case, our analysis shows that it is possible to read Heisenberg's paper as 
providing a complete (if limited) calculational scheme, the results of applying 
which are precisely those of standard quantum mechanics. Thus a second aim of 
our paper is to stress both the correctness and the practicality of what we con- 
jecture to be Heisenberg's calculational scheme, and to stimulate a re-appraisal 
of the possibility of including at least some discussion of it in undergraduate 
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courses. 



2 Heisenberg's 'transition amplitude' approach 
2.1 Quantum kinematics 

As is well known, Heisenberg begins his paper with a programmatic call to 'dis- 
card all hope of observing hitherto unobservable quantities, such as the position 
and period of the electron' 20 , and instead to 'try to establish a theoretical quan- 
tum mechanics, analogous to classical mechanics, but in which only relations 
between observable quantities occur'. As an example of such latter quantities, 
he immediately points to the energies W(n) of the (Bohr) stationary states, 
together with the associated (Einstein-Bohr) frequencies 21 

u)(n,n-a) = ~{W(n)-W{n-a)}, (1) 

noting that these frequencies, which characterise radiation emitted in the tran- 
sition n — > n — a, depend on two variables. An example of something he wishes 
to exclude from the new theory is the time-dependent position coordinate x(t). 
In considering what might replace it, he turns to the probabilities for transi- 
tions between stationary states. Consider a simple one-dimensional model of an 
atom consisting of an electron undergoing periodic motion, which is in fact the 
type of system studied by Heisenberg. For a state characterised by the label n, 
fundamental frequency u)(n) and coordinate x(n, t), one can represent x(n, t) as 
a Fourier series 

oo 

x(n, t) = \] A a (n)exp[iw(n)crf]. (2) 



=— oo 



According to classical theory, the energy emitted per unit time (that is, the 
power) in a transition corresponding to the ath harmonic uj(n)a is 22 



dE 



dt J n 37re e £ 



[co(n)a] A \X a (n)\ 2 . (3) 



In the quantum theory, however, the transition frequency corresponding to the 
classical l a>(n)a' is in general not a simple multiple of a fundamental frequency, 
but is given by thus u>(n)a is replaced by u>(n,n — a). Correspondingly, 
Heisenberg introduces the quantum analogue of X a (n), which he writes as 
X(n,n — a) 24 . Further, the left hand side of © has, in the quantum the- 
ory, to be replaced by the product of the transition probability per unit time, 
P(n, n — a), and the emitted energy huj(n, n — a); thus © becomes 

e 2 

P(n,n-a) = — — [ w (n,»- a)] 3 \X(n,n- a)\ 2 . (4) 

It is the transition amplitudes X(n,n — a) which Heisenberg fastens upon as 
being satisfactorily 'observable'; like the transition frequencies, they depend on 
two discrete variables 25 . 
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Equation (0} refers, however, to only one specific transition. For a full 
description of atomic dynamics (as then conceived), one will need to consider 
all the quantities X(n, n — a)exp[iaj(n, n — a)t]. In the classical case, the terms 
X a (n)exp[iLv(n)at] may be combined to yield x(t) via @. But in the quantum 
theory, Heisenberg says, a 'similar combination of the corresponding quantum- 
theoretical quantities seems to be impossible in a unique manner and therefore 
not meaningful, in view of the equal weight of the variables n and a [i.e. in the 
amplitude X(n,n — a) and frequency aj(n,n — a)].' 'However', he continues, 
'one may readily regard the ensemble of quantities X(n,n — a)exp[io;(n, n — a)t] 
as a representation of the quantity x(t) . . .' This is the first of Heisenberg's 
'magical jumps' - and certainly a very large one. Representing x(t) in this way 
seems to be the sense in which he considered himself to be 're-interpreting the 
kinematics'. 

Still concerned with the kinematics, Heisenberg immediately poses the ques- 
tion: 'how is the quantity x(t) 2 to be represented?' In classical theory, the 
answer is straightforward. From @ we obtain 

W)f = E E XaWX^n)^^'. (5) 

a 7 

Relabelling a + 7 as /?, 10 becomes 

[z(t)] 2 = EW eM " )/3 * ( 6 ) 



where 

Yp{n) =J2 X <*(n) x 0-<*{n). (7) 

a 

Thus classically [x(£)] 2 is represented (via a Fourier series) by the set of quan- 
tities Yp(n)exp[iu)(n)(3i\, the frequency uj(n)(3 being the simple combination 
[u>(n)a + uj(n){(3 — a)]. In quantum theory, the corresponding representative 
quantities must be written as Y(n, n — /3)exp[iw(n, n — f3)t], and the question is: 
what is the analogue of 10? 

The crucial difference in the quantum case is that the frequencies do not 
combine in the same way as the classical harmonics, but rather in accordance 
with the Ritz combination principle: 

u>(n, n — a) + u>(n — a,n — (3) — u)(n, n — f3), (8) 

which is of course consistent with 10. Thus in order to end up with the partic- 
ular frequency Lo(n,n — /3), it seems 'almost necessary' (in Heisenberg's words) 
to combine the quantum amplitudes in such a way as to ensure the frequency 
combination (JHJ); that is, as 

Y(n, n - p) e Mn,n-0)t = ^ x(n, n- a)e iuj{n ' n - a)t X(n -a,n- py^-^-P)* . 

a 

(9) 
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Cancelling the exponentials on both sides of © we are left with 



Y(n, n 



n,n — a)X(n — a, n — 0), 



(10) 



which is Hciscnberg's law for multiplying transition amplitudes together. 

He indicates the simple extension of the rule to higher powers [x(i)] n , but 
at once notices that a 'significant difficulty arises, however, if we consider two 
quantities x(t), y{t) and ask after their product x(t)y(t). . . . Whereas in classical 
theory x(t)y(t) is always equal to y(t)x(t), this is not necessarily the case in 
quantum theory'. This 'difficulty' clearly unsettled Heisenberg: but it very 
quickly became clear that the non-commutativity (in general) of kinematical 
quantities in quantum theory was the really essential new technical idea in the 
paper. Born recognised (|1U|) as matrix multiplication (something unknown to 
Heisenberg in July 1925), and he and Jordan rapidly produced Ref. 4, the first 
paper to state the fundamental commutation relation (in modern notation) 



Dirac's paper followed soon after 5 , and then the 'three-man' paper of Born, 
Heisenberg and Jordan 6 . 

The economy and force of Heisenberg's argument in reaching l|10|l is surely 
very remarkable - and it seems at least worth considering whether presenting it 
to undergraduates might not help them to understand the 'almost necessity' of 
non-commuting quantities in quantum theory. 

2.2 Quantum dynamics 

Having identified the transition amplitudes X(n, n — a) and frequencies to(n, n — 
a) as the 'observables' with which the new theory should deal, Heisenberg now 
turns his attention to how they may be determined 'from the given forces of the 
system' - that is, by the dynamics. In the Old Quantum Theory, he notes, this 
would be done in two stages: (a) by integration of the equation of motion 



and (b) by determining the constants of periodic motion through the 'quantum 
condition' 



where the integral is to be evaluated over one period of the motion. As regards 
<|12l) , Heisenberg says that it is 'very natural' to take over the classical equation 
of motion into quantum theory, replacing the classical quantities x{t) and f(x) 
by their kinematical re- interpretations 26 , as in the previous section (or, as we 
would say today, by taking matrix elements of the operator equation of motion 
x + f(x) =0). He notes that in the classical case a solution can be obtained 
by expressing x(t) as a Fourier series, insertion of which into the equation of 



xp — px = ih. 



(11) 



x + f(x) = 0, 



(12) 




(13) 
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motion leads (in simple special cases) to a set of recursion relations for the 
Fourier coefficients. In the quantum theory, Heisenberg says, 'we are at present 
forced to adopt this way of solving equation (|T2|) [his equation H(ll)]... .since 
it was not possible to define a quantum-theoretical function analogous to the 
[classical] function x{n,t).' 1 In section 3 we shall consider the simple example 
(the first of those chosen by Heisenberg) f{x) — ojqX + Ax 2 , obtaining the 
appropriate recursion relations in both the classical and the quantum cases. 

A 'quantum-theoretical re-interpretation' of (| 1 31) is now required, in terms 
of the transition amplitudes X(n,n — a). In the classical case, insertion of 
into ifHfl) gives 



mx 2 dt = 2ttto V \X a (n)\ 2 a 2 u(n) (14) 



a— — oo 



using X a (n) = [X- a (ri)]*', {T3J| is H(14). Heisenberg argues that does 
not sit well with the Correspondence Principle, since the latter should only 
determine J up to an additive constant (times h). Setting (|14|) equal to nh, he 
converts it to the form (H(15)) 

oo , 

/i = 2ttto a—{a\X a {n)\ 2 u{n)), (15) 

a— — oo 

which determines the X a (n)'s only to within a constant; the summation can 
alternatively be written as over positive values of a, replacing 27rm by Anm. 
In another crucial jump, Heisenberg now replaces the differential in (|15|) by a 
difference, giving 

oo 

h = 4nm ^~~^{| X (n + a, n)\ 2 u(n + a, n) — \X(n,n — a)\ 2 uj(n, n — a)}, (16) 

which is H(16) 27 . As he later recalled, he had noticed that 'if I wrote down this 
[presumably (|15|l above] and tried to translate it according to the scheme of 
dispersion theory, I got the Thomas-Kuhn sum rule [which is equation 116|l 28 ' 29 ]. 
And that is the point. Then I thought, "That is apparently how it is done" ' 30 . 

By 'the scheme of dispersion theory', Heisenberg is referring to what Jammer 31 
calls Born's correspondence rule, namely 32 

4-> *(n) - $(n - a), (17) 

on 

or rather to its iteration in the form 33 
d§(n, a) 



dn 



$(n + a,n) - <f>(n, n - a) (18) 



as used in the Kramers- Heisenberg theory of dispersion 34 ' 35 . It took Born only 
a few days to show that Heisenberg's quantum condition (|16|l was in fact the 
diagonal matrix element of (f 1 1 1> - and to guess 36 that the off-diagonal elements of 
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xp—px were zero, a result which was shown to be compatible with the equations 
of motion in Born and Jordan's paper 4 . 

At this point, when we have momentarily advanced a little beyond July 
1925, it may be appropriate to emphasize that Heisenberg's transition ampli- 
tude X(n, n — a) is indeed precisely the same as the quantum-mechanical matrix 
element (n — a|i|n), where \n) is the exact eigenstate with energy W(n). The 
relation of l|16fl to the fundamental commutator is briefly recalled in Ap- 
pendix A. 

Returning to the development of Heisenberg's paper, he summarizes the state 
of affairs reached so far by the statement that equations 1)1 2|l and 'if soluble, 
contain a complete determination not only of the frequencies and energies, but 
also of the quantum theoretical transition probabilities'. We draw attention to 
the strong claim here: that he has arrived at a new calculational scheme, which 
will completely determine the observable quantities. Let us now see in detail 
how this scheme works, for the case of a harmonic oscillator perturbed by an 
anharmonic force of the form Ax 2 per unit mass 37 . 

3 Heisenberg's calculational scheme, and its de- 
tailed working-out for the 'Ax 2 ' anharmonic 
oscillator 

3.1 Recursion relations in the quantum case 

The classical (and, by assumption, quantum-mechanical) equation of motion is 

x + uj 2 x + Xx 2 = 0. (19) 

Departing now from the order of Heisenberg's presentation, we shall begin by 
showing how - as he states - (I19fl leads to recursion relations for the transition 
amplitudes X(n,n — a). The (n, n — a) representative of the first two terms in 
1)19(1 is straightforward, being just 

[-Lu 2 (n, n-a)+ u 2 ]X(n, n - a)e iu{n ^- a)t , (20) 

while that of the third term is, by the rule (|10(l . 

A ^ X(n, n - (3)X(n - 0, n - a)e^ n - a ^ . (21) 

P 

The (n, n — a) representative 38 of (|19fl therefore yields 39 

[uj 2 - Lu 2 (n, n - a)]X(n, n - a) + A X(n, n - (3)X(n - /3, n - a) = 0, (22) 



which generates a recursion relation for each value of a (a = 0, ±1, ±2, . . .). For 
example, setting a — we obtain 

uj 2 X(n, n)+X[X(n, n)X(n, n)+X(n, n-l)X(n-l, n)+X(n, n+l)X(n+l, n)+. . .]. 

(23) 
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No general solution for this infinite set of non-linear algebraic equations seems 
to be possible, so - following Heisenberg - we turn to a perturbative approach. 



3.2 Perturbation theory 

To make the presentation self-contained, we need to include a certain number 
of ancillary results. Heisenberg begins by considering the perturbative solution 
of the classical equation (|12|) . He writes the solution in the form 

x(t) = Aao + ai cosojt + Xa 2 cos 2ujt + A 2 03 cos 2>u>t + . . . + \ a ~ 1 a a cos auit + . . . 

(24) 

where the coefficients a a , and w, are themselves to be expanded as a power 
series in A, the first terms of which are independent of A 38 : 

a a = 4 0) + a4 1} + A 2 a< 2) + . . . (25) 

o x = a^+Xa[ 1) +\ 2 a[ 2) + ... (26) 

and so on, and 

lu = uj a + Xlu (1) + AV 2) + . . . . (27) 

Straightforwardly inserting l|24|l into (|12ll . using standard trigonometric iden- 
tities, and equating to zero the terms which are constant, and which multiply 
coswt, cos2wt, etc, one obtains 

(constant) A{w 2 a + \^a\ + [A 2 (a^ + ^a 2 ) + ...]} = (28) 

(coswi) (-uj 2 +LLtl)a 1 + [\ 2 (a 1 a 2 + 2a a 1 ) + ...}=0 (29) 
(cos2wf) A{(-4w 2 +w 2 )a 2 + -a\ + [A 2 (aia 3 + 2a a 2 ) +...]} = (30) 

(cos3uit) \ 2 {(-9uj 2 + u;l)a 3 + a 1 a 2 + [\ 2 (a 1 a 4 + 2a a 3 ) + ...}} = (31) 

and so on, where the dots stand for higher powers of A. Dropping all the terms 
multiplying A 2 (and higher powers), - (|3*T)l become (for A 7^ and a\ 7^ 0) 

^a + ^a 2 = (32) 

i-^+uj 2 ) = (33) 

(-4j 2 + uj 2 )a 2 + \\ = (34) 

(-9ui 2 + w 2 )a 3 + a x a 2 = (35) 

which are the same as H(18) 41 . The 'lowest order in A' solution is now obtained 
from (|32|1 - i|35[) by setting u) = cjq, and replacing each a a by the corresponding 
one with a superscript l ( 0) ' (c.f. (J5SJ, ffify). 

Turning now to the quantum case, Heisenberg proposes to seek a solution 
analogous to Q24JI. Of course, it is now a matter of using the 'representation' 
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of x(t) in terms of the quantities {X(n, n — a)exp[icj(n, n — a)t]}. But it seems 
reasonable to assume that, as the index a increases away from zero, in integer 
steps, each successive amplitude will (to leading order in A) be suppressed by an 
additional power of A, as in the classical case. Thus Heisenberg suggests that, 
in the quantum case, x(t) should be represented by terms of the form 

Xa(n, n), a(n, n — 1) cos u)(n, n — l)t, Aa(n, n — 2) cos u>(n, n — 2)t, . . . , 
A Q_1 a(n, n — a) cos ui(n, n — a)t, . . . , (36) 

where, as in (|2SJ) - (|27j) . 

a(n, n) = a (0) (n, n) + Aa (1) (n, n) + A 2 a (2) (n, n) + . . . (37) 

a(n,n- 1) = a (0) ( n, n — 1 ) + Aa (1) (n, n - 1) + A 2 a (2) (n, n - 1) + . . . (38) 
and so on, and 

uj(n, n - a) = w (0) (n, n - a) + Aw (1) (n, n - a) + A 2 w (2) (n, n - a) + . . . (39) 

As Born and Jordan pointed out 4 , some use of 'correspondence' arguments has 
been made here, in assuming that, as A — > 0, only transitions between adjacent 
states are possible (we shall return to this point in section 3.3). 

Heisenberg now simply writes down what he asserts to be the quantum 
versions of (|32|l - Ij35(l . namely 42 

ula(n,n) + ^[a 2 (n + l,n) + a 2 {n,n- 1)] = (40) 

- 0J 2 {n,n - 1) +w 2 = (41) 
[-uj 2 (n,n- 2) + ujl]a(n,n - 2) + -[a(n,n- l)a(n -l,n- 2)] = (42) 

[— uj 2 (n 7 n — 3) + cjQ]a(n, — 3) + -[a(n, n — l)a(n — 1, n — 3) 
+a(n, n - 2)a(n - 2, n - 3)] = 0. (43) 

The question we now want to address is: how did Heisenberg arrive at equations 
(HJ - (E2J? 

Tomonaga 8 , having derived H22(l . then proceeds to discuss only the A — > 
limit - i.e. the simple harmonic oscillator, a special case to which we shall return. 
The only other authors, to our knowledge, who have offered a discussion of the 
presumed details of Heisenberg's calculations are 43 Mehra and Rechenberg 11 . 
They suggest that Heisenberg guessed how to 'translate', 're-interpret' or 're- 
formulate' (their words) the classical equations - (|35|) into the quantum 
ones (|4Uf) - (|43|l . in a way that was consistent with his multiplication law (|10|) . 
Although such 'inspired guesswork' was undoubtedly necessary in the stages 
leading up to Heisenberg's paper of July 1925, to us it seems more plausible 
that by the time of the paper's final formulation, Heisenberg realised that he 
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had a calculational scheme in which guesswork was no longer necessary, and in 
which (|i0"|) - in particular, could be derived. 

Unfortunately, we know of no documentary evidence which can directly prove 
(or disprove) this suggestion. But we think there is some internal evidence for it. 
In the passage to which attention was drawn earlier, Heisenberg firmly asserts 
that the formalism he has set up constitutes a complete scheme for calculating 
everything that needs to be calculated. It is hard to believe that Heisenberg did 
not realise that it led directly to equations lj4TJI) - |@3J|, without the need for any 
'translations' of the classical relations (of course, the latter were a nice check on 
the reasonableness of the quantum relations). At all events, it is the case that 
l|4UI) - <|43[) can be straightforwardly derived, as we shall now see. 

In order to apply the ansatz (|36|l to (|22|l . we need to relate the amplitudes 
X(n, n — a) to the corresponding quantities \ a ~ 1 a(n, n — a). We first note that 
in the classical case 

X a (n) = X*_ a (n) (44) 

since x(t) in © has to be real. Consider, without loss of generality, the case 
a > 0. Then the quantum-theoretical analogue of the left hand side of 144|) is 
X(n,n — a), while that of the right hand side is X*(n — a,n) (see Ref. 24). 
Hence the quantum-theoretical analogue of (|44|) is 

X(n,n — a) = X*(n — a, n), (45) 

which of course is nothing but the relation (n — a\x\n) — (n\x\n — a)* for the 
Hcrmitian observable x. Although X{n, n — a) can in principle be complex (and 
he twice discusses the significance of the phases of such amplitudes), Heisenberg 
seems to have assumed (as is certainly plausible) that in the context of the 
classical cosine expansion H24JI . and the corresponding quantum terms Q'lfijl . the 
X(n, n — a)'s should be chosen to be real, so that l|45|l becomes 

X(n, n — a) = X{n — a, n) (46) 

- that is, the matrix with elements {X(n,n — a)} is symmetric. Consider now 
a typical term of Ip?^ 

A Q ~ 1 a(n, n — a) cos[w(n, n — a)t] 
= — - — a(n, n — a){exp[ia;(n, n — a)t] + exp[— iu(n, n — ce)t]} — 

\ a — 1 

= — - — a(n, n — a){exp[iu;(ri, n — a)t] + exp[iw(n — a, n)t]} (47) 

using Lu(n,n — a) = —uj{n — a,n) from Q. Assuming that a(n,n — a) — 
a(n — a, n) as discussed for (|46f) . we see that it is consistent to write 

^a — 1 

X(n, n — a) = — - — a(n, n — a) (48) 

for positive values of a, and in general 

AH" 1 

X(n, n— a) — a(n,n — a), a^O. (49) 
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The particular case a = is clearly special, with X(n, n) = Xa(n, n). 

We may now write out the recurrence relations 1221) explicitly for a = 
0, 1,2,.. ., in terms of the a(n, n — a) rather than the X(n, n — a). We shall 
include terms up to and including A 2 . 

For a = we obtain 

\{u 2 a(n } n) + -{a 2 (n + 1, n) + a 2 (n, n — 1)] + 

X 2 [a 2 (n,n) + ^(a 2 (n + 2,n)+ a 2 (n,n- 2))]} = 0. (50) 

We note the connection with (|28|l . and that (|50|l reduces to Heisenberg's (|40|l 
when the A 2 term is dropped. Similarly, for a = 1 we obtain 



(— uj 2 (ii, n — 1) + a;g)a(n, n — 1) + \ 2 {a(n, n)a(n, n — 1) + a(n, n — l)a(n — 1, n — '. 
— [a(n, n + l)a(n 

(c.f. P^ll): for a = 2 



— [a(n, n + l)a(n + 1, n — 1) + a(?i, n — 2)a(rt — 2, n — 1)]} = 



A{(— u! 2 (n, n — 2) + aj 2 )a(n, n — 2) + — a(n, n — l)a(n — 1, n — 2) + 

A 2 [a(n, n)a(n, n — 2) + a(n, n — 2)a(n — 2, rt — 2) + — a(ri, n + l)a(ri + 1, rt — 2) 

+ ia(n,n- 3)a(n- 3,n- 2)]} = (52) 

(c.f. (J201)); and for a = 3 (c.f. l|5T]l) 

A 2 {(— cj 2 (n, n — 3) + cjQ)a(n, n — 3) + 

— [a(n, n — l)a(n — 1, n — 3) + a(«, w — 2)a(n — 2, n — 3)] + 

A 2 [a(n, n)a(n, n — 3) + a(n, n — 3)a(n — 3, n — 3) + -a(«, n + l)a(n + 1, rt — 3) + 
-a(n, n - A)a{n - 4, n - 3)]} = 0. (53) 

On dropping the terms multiplied by A 2 , (|50|l - (|53|l reduce to Heisenberg's <|4(JI) 
- (|43|l . This appears to be the first published derivation of these equations. 

In addition to these recurrence relations which follow from the equations of 
motion, we also need the perturbative version of the quantum condition (|16|) 44 . 
Including terms with a A 2 power, consistent with l|5U|) - (|53|l . (|16fl becomes 

= a 2 (n + 1, n)uj(n + 1, n) — a 2 (n, n — l)cj(n, n — 1) 

Tim 

+X 2 [a 2 (n + 2, n)uj(n + 2, n) - a 2 (n, n - 2)u(n, n - 2)]. (54) 
We are now ready to obtain the solutions. 
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3.3 The lowest-order solutions for the amplitudes and fre- 
quencies 

We begin by considering the lowest-order solutions, in which all A 2 terms are 
dropped from equations H5U|) - (|54|l . and all quantities (a's and w's) are replaced 
by the corresponding ones with a superscript ^ (c.f. I|37|) - l|39fl 45 . In this case, 
ifST]) reduces to 

[-(J°\n, n - l)) 2 + w 2 ]a (0) (n, n - 1) = 0, (55) 
so that assuming a^(n, n — 1) 7^ we obtain 

LO (0 \n,n- 1) = uj q (56) 
for all n. Substituting (|56(l into the lowest-order version of (|54() we find 

— — = [ a f) (n + 1, 7^] 2 - [a (0) (n, n - l)] 2 . (57) 

The solution of this difference equation is 

[a {0) (n 7 n - l)] 2 = — — (n + constant) , (58) 
Trmujo 

as given in H(20) 45 . To determine the value of the constant, Heisenberg used 
the idea that in the ground state there can be no transition to a lower state. 
Thus 

[a (0) (0,-l)] 2 = (59) 

and the constant is determined to be zero. This then gives (up to a convention 
as to sign) 

a (0 \n,n- 1) = (60) 

where 

P = {h/nmLUo) 1 ^. (61) 

Equations (|56|l and H6()J) are Heisenberg's first results, and they are in fact 
appropriate to the simple (unperturbed) oscillator. We can check (|60|) against 
standard quantum mechanics via 

a [0) {n,n- 1) = 2X^(n,n-l) = 2 (n~l\x\n) Q (62) 

where the states |n)o are unperturbed oscillator eigenstates. It is well known 
that 46 



o(ti - l\x\n) = y 2muJa J ^ ( 63 ) 
which agrees with l|6U|) . using l|61|) . A similar treatment of i|5U[) leads to 

a^{n,n) = —^{2n + l). (64) 
4w 
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Turning next to equation l|52|) . the lowest-order form is 



(-[J Q Hn,n-2)] 2 +Lol)a^\n,n-2) + ^a {0) (n,ri-l)a^ \n-l,n~2) = 0. (65) 

Now the composition law must be true for the lowest-order frequencies, so we 
have 

uj {0) {n,n~2) = u {a) (n,n-l)+uj ( - 0) {n-l,n-2) 

= 2uj (66) 

using (|56|l . and in general 

uj {0) {n,n- a) = auj (a = 1, 2, 3, . . .)■ (67) 
Using (jnOJ, © and we then obtain 





a (0) (n,n-2) = ±Wn(n - 1). (68) 



A similar teatment of l|53|) yields 



/3 : 
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a<°>(n,n - 3) = ^ y/n(n - l)(n - 2). (69) 
Consideration of the general lowest-order term in (|22ll shows that 



^Hn,n-a) = A a -^ ] j ¥ ^ : (70) 

where A a is a numerical factor depending on a; (|70|l is equivalent to H(21). 

At this stage it is instructive to comment briefly on the relation of the 
above results to those which would be obtained in standard quantum-mechanical 
perturbation theory. At first sight, it is surprising to see non-zero amplitudes 
for two-quantum ( 1)68(1 ). three-quantum 1 1)69(1 ) or a-quantum (((ZOJ) transitions 
appearing at 'lowest order'. But we have to remember that in Heisenberg's 
perturbative ansatz (|36|l . the a-quantum amplitude appears multiplied by a 
factor A Q_1 . Thus, for example, the 'lowest order' two-quantum amplitude is 
really Aa^°^(n,7i — 2), not just a^°\n,n — 2). Indeed, such a transition is to 
be expected precisely at order A 1 , in conventional perturbation theory. The 
amplitude is (n — 2\x\n) where, to order A, 

i \ i \ , 1 \ o(k\x 3 \n} 

|n) = k)o + -mAg lfc) . (71) 

The operator x 3 connects |rt)o to \n + 3)o, \n + l)o, \n — l)o, and \n — 3)o, and 
similar connections occur for q(ti — 2|, so that a non-zero O(A) amplitude is 
generated in (n — 2|i|n). 
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It is straightforward to check that is indeed correct quantum-mechanically; 
but it is more tedious to check (|69|l . and distinctly unpromising to contemplate 
checking l|7U|) by doing a conventional perturbation calculation to order a — 1. 
For this particular problem, the 'improved' perturbation theory represented by 
OSl) is clearly very useful. 

Having calculated the amplitudes for this problem to lowest order, Heisen- 
berg next considers the energy. Unfortunately he again gives no details of his 
calculation, beyond saying that he uses the classical expression for the energy, 
namely 

W = \^nx 2 + ^itlujqX 2 + ^m\x 3 . (72) 

It seems a reasonable conjecture, however, that he replaced each term in <|72ll 
by its 'ensemble of representatives'. Thus x 2 , for example, is replaced by the 
ensemble of terms of the form 

X(n, n - f3)X(n -p,n- a ) e ^(n,n-a)t ^ (73) 



according to his composition law (|1C)|> . Similarly, x 2 is replaced by terms of the 
form 

k>(n, n - f3)X(n, n - fif&^-Miufa - /3, n 



= u)(n, n — (3)uj{n — a,n — (3)X(n, n — f3)X(n 


using to(n, m) = —Lj(m, n). In general then, in Hcisenberg's scheme, the 'energy 
representatives' will be of the form 

W(n,n-a)e iuJ< - n ' n - a)f . (75) 

It follows that if energy is to be conserved (i.e. time-independent) the off-diagonal 
elements must vanish: 

W(n, n - a) = 0, a ^ 0. (76) 

The term a = is time- independent, and may be taken to be the energy in 
the state n. The crucial importance of checking the condition Il76|l was clearly 
appreciated by Heisenberg. 

To lowest order in A, the last term in (|72(l may be dropped. Furthermore, 
referring to 13(i|) . the only A-independent terms in the X-amplitudes are those 
involving one-quantum jumps such as n — ► n — 1, corresponding in lowest order 
to amplitudes such as X^(n, n — 1) = ^a^(n, n — 1). This means (referring 
to O and (J71J) that the elements W(n, n), W(n, n-2) and W(n, n + 2) , and 
only these W(n, m)'s, are independent of A when evaluated to lowest order. In 
Appendix B we show that W(n, n — 2) vanishes to lowest order, and W(n, n + 2) 
vanishes similarly. Thus, to lowest order in A, the energy is indeed conserved (as 



- a)X{n- (3,n- ay^-P 
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Heisenberg notes), and is given (using (|73|l and 174|) with a = and (3 = ±1) 

by 

W(n,n) = ^m[uj^{n,n-\)] 2 [X^{n,n-\)] 2 + 

n + 1, n)] 

+ imc 2 [X(°)(n,n - l)] 2 + ^[A^n + l,n)] 2 

= (n+~)7kJ (77) 

using 1)56(1 . and This is the result given by Heisenberg in H(23). 

These 'lowest order' results are the only ones Heisenberg reports for the Ax 2 
force. We do not know whether he carried out higher-order calculations for this 
case or not. What he writes, next, is that the 'more precise calculation, taking 
into account higher order approximations in W, a and uj will now be carried 
out for the simpler example of an anharmonic oscillator x + uj 2 + Ax 3 = . ' 
This case is slightly simpler because only the 'odd' terms are present in l|36[) 
(i.e. ai,Aa 3 ,A 2 a 5 , etc). The results Heisenberg states for the 'Ax 3 ' problem 
include terms up to order A in the amplitudes, and terms up to order A 2 in the 
frequency uj{n,n — 1), and in the energy W. Once again, he gives no details 
of how he has done the calculations. We believe there can be little doubt that 
he went through the algebra of solving the appropriate recurrence relations, up 
to order A 2 in the requisite quantities. As far as we know, the details of such 
a calculation have not been given before, and we therefore feel that it is worth 
giving them here, as being of both pedagogical and historical interest. In the 
following section we shall obtain the solutions for the Ax 2 force (up to order A 2 ) 
which we have been considering hitherto, rather than start afresh with the Ax 3 
one. The procedure is of course the same for both. 

Before leaving the 'lowest order' calculations, we address a question which 
may have occurred to the reader: given that, at this stage in his paper, the 
main results actually relate to the simple harmonic oscillator rather than to the 
anharmonic one, why did Heisenberg not begin his discussion of 'toy models' 
with the simplest one of all, namely the simple harmonic oscillator? And indeed, 
is it not possible to apply his procedure to the SHO, without going through the 
apparent device of introducing a perturbation, and then retaining only those 
parts of the solution which survive as the perturbation vanishes? 

For the unperturbed oscillator (the SHO), the equation of motion is of course 
x + u)qX = 0, which yields 

[luI - uj 2 (n,n- a)]X(n,n- a) = (78) 

for the amplitudes X and frequencies to. It is also reasonable to retain the 
'quantum condition' lltjl) . since this is supposed to be true whatever the partic- 
ular force law. If we assume that the only non-vanishing amplitudes are those 
involving adjacent states (because, for example, in the classical case only a sin- 
gle harmonic is present 47 ), then - remembering that X(n, n — 1) = ha{n, n — 1) 
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- (tTol) and l(7Sll reduce to l(57|l and respectively, and we quickly recover 
the previous results. Thus we see that this is indeed an efficient way to solve 
the quantum SHO 48 . For completeness, however, it would be nice not to have 
to make the 'adjacent states' assumption; Born and Jordan 4 showed how this 
could be done, but their argument is somewhat involved. Soon thereafter, of 
course, the wave-mechanics of Schrodinger, and the operator approach of Dirac, 
provided the derivations used ever since. 

We now turn to the higher order corrections, for the Aa; 2 force. 



3.4 The solutions up to and including A 2 terms 

Consider first equation l|51|l . retaining terms of order A but no higher powers. 
We set 

tu(n,n-l) = u>o + \cu (1) (n,n-l) (79) 
a{n,n-l) = a (0) (n, n - 1) + Aa (1) (n, n - 1) (80) 

and find 

2\LU uj {1) (n,n- l)a (0) (n,n- 1) = 0, (81) 

so that 

uj {1 \n,n~ 1) = 0. (82) 

Consideration of equation H51J) up to terms of order A 2 , employing equations 
<|6U|) . i|64l) and (|68[1 for the zeroth order amplitudes, gives the 0(A 2 ) correction 
to u(n,n- 1) (c.f. (E3)): 

0J {2 \n,n-l) = -^n. (83) 

The corresponding corrections to a(n, n — 1) are found from the quantum 
condition l|16|) . To order A we set 

a{n + l,n) = a (0) (n + l,n) + Xa (1) (n+ l,n) (84) 

as in H8U[I. and find 

V^TTa (1) (n + 1, n) - V™a (1) (n, n - 1) = 0. (85) 

This equation has the solution a^{n, n— 1) = constant /y/n, but the condition 
a^(0, —1) = (c.f. (|59|) ~) implies that the constant must be zero, and so 

a {1) {n,n- 1) = 0. (86) 

In a similar way, to order A 2 we obtain 

V^+la {2) (n + 1, n) - V^a (2) (n, n - 1) = i^(2n + 1), (87) 
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which has the solution 

a ( - 2 \n,n-l) = ^ 1 n^. (88) 



n 



We can now find the higher order corrections to a(n, n) by considering equa- 
tion JSUJ). We obtain 



B 4 

a [L) (n,n) = 0, a w (n,n) = -^-g (30n 2 + 30n + 11). 



Similarly, we find from (|5 



a (1) (n,n-2) = 0, a (2) (n,n-2) = ^(2n - l)^n{n - 1), (90) 

oZLOq 

where we have used 

w (2) (n,n-2) = w (2) (n,n-l)+w (2) (n-l,n-2) 

° '' (2n-l). (91) 
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These results suffice for our purpose. When n is large, they agree with those 
obtained for the classical 'Ax 2 ' anharmonic oscillator using the method of suc- 
cessive approximations 49 . As an indirect check of their quantum-mechanical 
validity, we now turn to the energy, evaluated to order A 2 . 

Consider first the (n,n) element of ^muigX 2 . This is given, to order A 2 , by 

imw 2 {i[(a(°)(n,n- l)) 2 + (o<°>(n,n+ l)) 2 ] + y [4(a<°> (n, n)f + 

2a (2) (n,n- l)a (0) (n - l,n) + 2a (2) (n,n+ l)a (0) ( n + 1, n) + 
( a (°) (n, n - 2)) 2 + (a(°) (n, n + 2)) 2 ]} 

- \™^{^{n +h + +n+ 11/30)}. (92) 



Similarly, using l|74|) up to order A 2 , with a = 0, the (n, n) element of ^mx is 



found to be 



1 2 (3 2 1 5/3 4 A 2 
- roWo { y (n+-)- _ 



-m^{^(n + -) - ^(n 2 + n + 11/30)}. (93) 



Finally we need to consider the (n,n) element of the potential energy ^mXx 3 . 
To obtain the result to order A 2 , we need compute the (n, n) element of x 3 only 
to order A. Using 



x 3 



(n, n) = ^2 n ~ a ) x ( n -ot,n- 8)X(n - 3, n) (94) 

a p 
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one finds that there are no zeroth order terms, but twelve terms of order A 
(recall that amplitudes such as X(n,n) and X(n,n — 2) each carry one power 
of A) . Evaluating these terms using the previously obtained results gives 

~~24^r ( n +" + 11 /30) (95) 

for this term in the energy. Combining 192J1 . (|93[1 and H95[l then gives, for the 
energy up to order A 2 , 

1 5X 2 h 2 

W(n,n) = {n+~)%u -— j(n 2 + n + 11/30), (96) 

£ IZTYILJq 

a result 50 which agrees with classical perturbation theory when n is large 51 , 
and is in exact agreement with standard second-order perturbation theory in 
quantum mechanics 52 . 

As mentioned earlier, Heisenberg does not give results for the 'Ax 2 ' force 
beyond zeroth order; he does, however, give the results for the 'Aa; 3 ' force up 
to and including A 2 terms in the energy, and A terms in the amplitudes. By 
'the energy' we mean, as usual, the (n, n) element of the energy operator, which 
as noted in section 3.3 is independent of time. One should also check that the 
off-diagonal elements W(n,n — a) vanish (see equation (|76(l ^l. These are the 
terms which would (if non-zero) carry a periodic time-dependence, and in his 
paper Heisenberg says that 'I could not prove in general that all periodic terms 
actually vanish, but this was the case for all the terms evaluated'. We do not 
know how many off-diagonal terms W(n, n — a) he did evaluate, but he clearly 
regarded their vanishing as a crucial test of the formalism. In Appendix B we 
outline the calculation of all off-diagonal terms for the Ax 2 force, up to order A, 
as an example of the kind of calculation Heisenberg probably did, finishing it 
late one night on Heligoland 53 . 



4 Conclusion 

We have tried to remove some of the barriers to understanding Heisenberg's 
1925 paper, by providing (apparently for the first time) details of calculations 
of the type we believe he probably carried out. We hope that more people will 
thereby be encouraged to appreciate this remarkable paper. 

The fact is, Heisenberg's 'amplitude calculus' works: at least for the simple 
one-dimensional problems on which he tried it out, it is an eminently practical 
procedure, requiring no sophisticated mathematical knowledge to implement. 
Since it uses the correct equations of motion, and incorporates the fundamen- 
tal commutator l|llf) via the 'quantum condition' (|16[) . the answers obtained are 
completely correct, in the sense of agreeing with conventional quantum mechan- 
ics. 

We believe that Heisenberg's approach, as applied to these simple dynamical 
systems, has considerable pedagogical value, and could usefully be included in 
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undergraduate courses on quantum mechanics. The multiplication law IjlUfl 
has a convincing physical rationale, even for those who (like Hcisenberg) do not 
recognize it as matrix muliplication; indeed, this piece of quantum physics could 
provide an exciting application for those learning about matrices in a concurrent 
mathematics course. The simple examples of <|1(J[) . in formulae such as <|22[1 
or the analogous one for the Ai 3 force, introduce the student directly to the 
fundamental quantum idea that a transition from one state to the other occurs 
via all possible intermediate states, something which can take time to emerge 
in the traditional wave-mechanical approach. The solution of the quantum 
SHO, sketched at the end of the previous section, is simple and elementary, in 
comparison with the standard methods. Finally, the type of perturbation theory 
employed here provides an instructive introduction to the technique, being more 
easily related to the classical analysis than is conventional quantum-mechanical 
perturbation theory (which students tend to find very formal) . 

It is of course true that many important problems in quantum mechanics are 
much more conveniently handled in the wave-mechanical formalism: unbound 
problems are an obvious example, but even the Coulomb problem required a 
famous tour de force from Pauli 54 . Nevertheless, a useful seed may be sown, so 
that when students meet problems involving a finite number of discrete states - 
for example, in the treatment of spin - the introduction of matrices will come as 
less of a shock. And they may enjoy the realisation that the somewhat myste- 
riously named 'matrix elements' of wave-mechanics are indeed the elements of 
Heisenberg's matrices. 



Appendix A: The quantum condition 
(1161) and xp — px = \h. 

Consider the (n,n) element of (xx — xx). This is 

X(n, n — a)iu>(n — a, n)X{n — a, n) 

a 

— \uj(n, n — a)X(n, n — a)X(n — a,n). (97) 

a 

In the first term of i|9 7p. the sum over a > may be re- written as 

-i^2uj{n,n~a)\X{n,n-a)\ 2 (98) 

using Lj(n, n — a) = —u>(n — a,n) from and X(n — a, n) = X*(n, n — a) from 
l|45(l . while the sum over a < becomes, similarly, 

i^2u}{n + a,n)\X(n + a, n)\ 2 (99) 



Q>0 
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on changing a to —a. Similar steps in the second term of (|97|) lead to the result 
(xx — xx)(n, n) — 2i [oj(n + a, n)\X(n + a, n)\ 2 

a>0 

~u>(n 7 n — a)\X(n, n — a)| 2 ] 
= 2i/i/(4?rm), (100) 

where the last step follows from the 'quantum condition' i|16[) . Setting p = mx 
we find 

(xp — px)(n, n) — ifi (101) 

for all values of n. This is the result which Born found 36 shortly after reading 
Heisenberg's paper. In the further development of the theory the value of the 
'fundamental commutator' xp — px, namely ih times the unit matrix, was taken 
to be a basic postulate. The sum rule (|16|) is then derived by taking the (n, n) 
matrix element of the relation [x, [H, x}] — h 2 /m. 



Appendix B: Calculation of the off-diagonal 
matrix elements of the energy W(n, n — a) 
for the Xx 2 force, up to order A. 

We shall show that, for a ^ 0, all the elements (n,n — a) of the energy 
operator Imi + Imwgi 2 + ^Xmx 3 vanish, up to order A. 

We begin by noting the qualitative point that, at any given order in A, 
only a limited number of elements W(n, n — 1), W(n, n — 2), . . . will contribute, 
since the amplitudes X(n, n — a) are suppressed by increasing powers of A 
as a increases. In fact, for a > 2 the leading power of A in W(n,n — a) is 
A"~ 2 , arising from terms such as X(n,n — l)X(n — l,n— a) and XX(n,n — 
l)X(n — l,n — 2)X(n — 2,n — a). Thus to order A we need only calculate 
W(n, n - 1), W(n, n - 2), W(n, n - 3). 

(a) W(n,n-l) 

There are four O(A) contributions to the (n, n — 1) element of \mio^x 1 : 

-mw 2 A{a (0) (n,n)a (0) (n,n- 1) + a (0) (n,n- l)a (0) ( n — 1, n — 1) + 

- [ fl (°) (n, n + l)a^ (n + 1, n - 1) + a (0) (n, n - 2)a^ (n - 2, n - 1)]} 
5 — 

= -— m\f3 3 ny/^. (102) 
There are two O(A) contributions to the (n,n — 1) element of ^mx : 

-- Am{w (0) (n, n + l)w (0) (n + 1, n - l)a (0) (n, n + l)a (0) (n + 1, n - 1) + 
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w (0) (n, n - 2)w (0) (n - 2, n - l)a (0) (n, n - 2)a (0) (n - 2, n - 1)} 

= -^m\(3 3 n^. (103) 

There are three 0(A) contributions to the (n, n — 1) element of |mAi 3 : 

— mA{a (n, n — 1) + 

a (0) (n, n- l)a (0) (n- l,n-2)a (0) ( n — 2, n — 1) 
W 0) (n, n + l)a(°'(n + 1, n)a (0) (n, n - 1)} 

= -m\(3 3 n^. (104) 
8 

The sum of (I102|) . (|1L)3|) and (|104|l vanishes, as required. 

(b) W(n,n-2) 

The leading contribution is independent of A. From the term ^moJ^x 2 it is 

-mu 2 a (0) (n, n - l)a (0) (n - 1, n - 2), (105) 
8 

which is cancelled by the corresponding term from ^mx . The next terms are 
0(A 2 ), for example from the leading term in the (n, n — 2) element of ^Xmx 3 . 

(c) W(n,n-3) 

There are two O(A) contributions from huiujQX 2 '. 



1 

1 



mw 2 A{a (0) (n, n - l)a (0) (n - 1, n - 3) + a (0) (n, n - 2)a (0) (rc - 2, n - 3)} 



m\(3 3 ^n(n -l)(n- 2). (106) 
There are two O(A) contributions from : 
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mA{w (0) (n, n - l)a (0) (n, n ~ l)w (0) (n - 1, n - 3)a (0) (n - 1, n - 3) 



W 0) (n, n - 2)a (0) (n, n - 2 V (0) (n - 2, » - 3)a (0) (n - 2, n - 3)} (107) 
= -— Xm/3 3 ^n(n - l)(n - 2). (108) 

There is only one O(A) contribution from ^mAa; 3 : 

-^■mAa (0) (n, n - l)a (0) (n - 1, n - 2)a (0) (n - 2, n - 3) 
= ^Xmf3 3 ^n(n- l)(n-2). (109) 
The sum of (I1U6|) . (|108|) and (|109|l vanishes, as required. 
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